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It is shown that the order parameter A induced in the normal part of superconductor-normal- 
superconductor proximity system is modulated in the magnetic field differently from vortices in bulk 
superconductors. Whereas A turns zero at vortex centers, the magnetic structure of these vortices 
differs from that of Abrikosov’s. 


The question of superconductivity induced in the nor¬ 
mal part (N) of superconductor-normal-superconductor 
(SNS) proximity system has recently been revived by ob¬ 
servations of vortices in N [1]. The order parameter A 
induced in N is not uniform even in zero field and is 
strongly suppressed nearly everywhere in N except the 
immediate vicinity of interfaces. Hence, the formal prob¬ 
lem of the order parameter distribution within vortices 
in N is qualitatively different from that of bulk supercon¬ 
ductors and so do the physical properties of “N-vortices”. 
These properties are of interest both for the basic physics 
and for applications, enough to mention wires in super¬ 
conducting magnets which are in fact SNS systems. 

Describing proximity effects, one encounters the ques¬ 
tion of the length scale on which the induced order pa¬ 
rameter varies. This problem is discussed in the first part 
of this paper for any field and temperature. In the fol¬ 
lowing part, a linear combination of the eigenfunctions 
of an equation for A is constructed to represent vortices 
in N. In fact, the seminal work of Abrikosov on type-II 
superconductors suggests the form of this combination 
[2]. The difference, though, is that Abrikosov combined 
eigenfunctions of the 1st Landau level, whereas in the 
problem of interest here these functions are different and 
more general. 

As mentioned, the induced A is strongly suppressed 
nearly everywhere in N except the vicinity of interfaces. 
Out of this vicinity, equations of superconductivity in N 
can be linearized. Formally, the situation is similar to 
that at the upper critical field iLc 2 5 where the magnetic 
field is uniform and the small A satisfies a linear equation 

-^2n2A = A, or U‘^A = k‘^A. (1) 


at any temperature T [3]. Here, H = V + 27TiA/(pQ, A 
is the vector potential, 0o is the flux quantum, and k‘^ = 
— 1/^^. Notwithstanding the form, this equation differs 
from the linearized Ginzburg-Landau equation (GL), in 
the latter the coherence length ^ diverges as T Tc. At 
Hc 2 and T ^ Tc, is finite and is found by solving 

the self-consistency equation of the theory 
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/3 = 1 + 2cjr . (2) 


Here, huu = 7rT(2n + 1) are Matsubara energies and r is 
the scattering time for non-magnetic impurities. Accord¬ 


ing to Helfand and Werthamer [3], 



where i = vr is the mean-free path. 

Eq.0 is equivalent the Schrodinger equation for a 
charge in uniform magnetic field; Hc 2 = cor¬ 

responds to the minimum eigenvalue. The correspond¬ 
ing eigenfunctions belong to the first Landau level. A 
linear combination of these functions, constructed by 
Abrikosov, represents the lattice of vortices [2]. 

The normal metal within the proximity system may 
have its own < T and Hc 2 ,n{T). We are interested 
here in the part of the phase diagram outside of the re¬ 
gion under Hc 2 ,n{T) (within this region the N part is 
superconducting and the proximity system should rather 
be called SS'S). In this broad domain, the induced super¬ 
conductivity is still described by Eqs. Q and ([^, how¬ 
ever with a more general S{H,T,r) [11 


S{H,T,r) = 
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Here erfcz = 2 e Solving the self- 

consistency Eq. (§ with the new S, one can evaluate 
^(iL, T, r) in any place of the {H,T) phase diagram. 

At Hc 2 ,n{T), = 0o/27riL i.e. k^/q^ = -1, and 

the parameter a = —1. Therefore, the self-consistency 
equation (§ in dimensionless form, 

-Inf-^ ^ ^(2n + 1) - A5 / ’ 

n=0 ^ \ / / 


(A = h/27rTcN^ is the non-magnetic scattering param¬ 
eter) should give Hc 2 ,n{T) if one sets a = —1 in S' of 
Eq. Solving this numerically for the clean limit one 
obtains the lower curve of Eig.[^ see Appendix A. 

If iL -> 0, the parameter a diverges, whereas u ^ 0. It 
is readily shown [5] that S of this case has a closed form: 

R kjf 

5(0,T,r) = gtanh-i-. (7) 

Solving numerically Eq. ([^ with S(0,T, r) one obtains 
the decay length ^ = 1/k of the order parameter in the 
normal part of proximity systems mm in zero field. 
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FIG. 1. (Color online) The lower curve is Hc 2 of the clean 
limit in units at this curve = —q^ and 

(7 = — 1. At the upper curve = 0 and a = —1/2. Between 
the upper and lower curves, k^ is negative. Above and to the 
right of the upper curve, k^ is positive. 


Thus, k‘^ = —< 0 at the curve Hc 2 whereas it must 
be positive in zero field at t > 1, where it describes A 
attenuation in the N phase. This suggests that a curve 
exists on the plane (H^T) where = 0 [4]. This ques¬ 
tion is addressed by setting a = —1/2 in S of Eq. @ 
and solving the latter for This curve evaluated 

numerically for the clean limit is the upper one in Fig.[^ 
The question then arises about behavior of the normal 
metal in the part of the (i7, T) plane where k‘^ is negative 
(between the curves of Fig.[^ and out of it where k^ > 0. 
To address this we look at eigenfunctions A{x,y) of the 
equation 11^A = k‘^A. Choosing A = —Hyx we have 


d 2 itH 




dy‘^ 


(8) 


The equation does not contain x explicitly, so that 


A = Aoe*P"x(2/) (9) 


with x{y) satisfying 

d'^x/dy'^-<f{y + p/<ffx = k'^x- ( 10 ) 

In terms of y = y p/q^ the general solution is: 


X = CiXi + C2X2 , 

Xi = 2(<^+i)/2e-9^s/V2 ^ 

X2 = ^ 

with arbitrary constants Ci ,2 and cr of Eq. The Her- 
mite functions H(cr, w) can be expressed in terms of the 
parabolic cylinder functions and reduce to Hermite poly¬ 
nomials for cr = 0,1, 2, ... [8]. 


X\ 



FIG. 2. (Color online) Xi(cr,y) for g = 1 in units ^'KTcjhv. 
Ordering curves by their left edges clockwise: g — (i4c2), 

—0.9 (with a minimum at the left), —0.5 (/c = 0), and 1. 


Note that x\ with cr being a negative integer are the 
harmonic-oscillator wave-functions that go to 0 as ^ ^ 
Too; these are the eigenfunctions of the Landau levels. 
We are interested here in the part of the phase diagram 
where cr > —1 where Xi(^) is real, diverges as ^ ^ —oc, 
and goes to 0 as ^ ^ +oc, see Fig.[^ For symmetric SNS 
systems, xi should be discarded. 

On the other hand, X 2 in general has both real and 
imaginary parts. An example of X 2 (—0.7, 1., ^) is shown 
in Fig.[^ Both real and imaginary parts diverge at large 
y. X 2 should be taken into account in finite samples, 
unlike the infinite ones where it should be disregarded. 







FIG. 3. (Color online) The function X 2 (cr, y) for cr = —0.7 
and q—\. Re[x 2 (y)] is even, Im[x 2 (y)] is odd function of y. 

Consider now a normal metal layer between two thick 
superconducting banks forming the SNS proximity junc¬ 
tion. The N slab of a thickness W is infinite in x and z 
directions whereas —W/2 < y < W/2. The temperature 
of the system Tcn < T < In zero field Eq. (§ de¬ 
scribes exponential decay of the order parameter away of 
interfaces, A ex cosh^/^. In the magnetic field along ^ 
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which is less than Hd^s^ the field is confined to the N do¬ 
main, whereas the S banks are in the Meissner state (as¬ 
suming the London penetration depth <C W). Since 
superconductivity is induced in N by proximity with S', 
one expects vortices to be nucleated within the N layer. 
In small enough fields, vortices should form a periodic 
chain in the slab middle at ^ = 0. 

The N slab is uniform in the x direction, so that the 
parameter p in the solution can take any value. Con¬ 
sider a linear combination 


A = e^^^X2{y+p/q^) + e-'*>"X2(y -p/g"), (12) 


where the overall constant Aq is omitted. It is clear that 
if A(xo,0) = 0, the zero should be repeated with the 
period Sxq = n/p. If the penetration depth into the S 
banks is small relative to IT, the flux quantization would 
have given SxqWH = 0o and the parameter p would be 


p = 


nWH 

00 


q^W 

2 


(13) 


Unlike the problem of Hc 2 ,n{T) where = hv/27TTcN 
was adopted as a natural unit length, it is convenient 
here to normalize lengths to IU/2, the half-width of the 
N layer. Then, the dimensionless p = q^. Since the RHS 
of Eq. (12) is dimensionless, we keep the same notation 
x^y^p^q as for their dimensional counterparts. 

The structure of the solution is illustrated in 

Fig.|^where the modulus \A{x^y)\ is plotted for IT = 2, 
cr = 0.2, g^=p = 2. As predicted, the distance be¬ 
tween singularities (vortices) is 6xo = n/p ^ 1.57. This 
means that the flux quantization holds indeed for vor¬ 
tices in the N layer, which is not self-evident in advance. 
The solution shown is normalized as to have A 1 at 
the interfaces y = ±1, i.e., Aq is set equal to inverse of 
the RHS of Eq. taken at ^ = 1. The phase of this 



FIG. 4. (Color online) Lin es o f the constant order parameter 
modulus according to Eq. (12) with a = 0.2, ^ = a/ 2, p = 2 


for —2 < X < 2 and —1 < y < 1. 


solution near one of the vortices is shown in Eig.[^ 

It should be noted that the form (12) of the order pa¬ 


rameter is not the only possibility. One can consider 
various linear combinations with different complex coef¬ 
ficients which all satisfy H^A = A. The choice of these 
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FIG. 5. (Color online) Contours of the constant phase with 
the step 7r/4 in clock-wise order for the vortex at ir ~ 0.6 and 
p = 0 of Fig.|^ The phase jumps by 27r at the straight ray 
from the vortex center to the left. 


coefficients is dictated by the boundary conditions (the 
form of the S banks and the distribution of the order 
parameter on these banks). An example of |A| for 

A = e^^-X2{y + p/q^) + 5(1 + i)e-^^^X2{y - p/q^) (14) 
is shown in Fig.In this example A is not normalized 



FIG. 6. (Color online) Contours of |A| = constant for the 
order parameter of Eq. (14) with a — —0.8 and p — q^ — 2. 


and it clearly does not satisfy the boundary conditions 
AN{y = il) = A^, but it shows qualitatively that this 
type of linear combinations might be useful in describing 
the proximity effect in asymmetric SNS' systems with 
different S-banks for which vortices in N tend to be closer 
to the bank with smaller order parameter. The value 
a = —0.8 is chosen to illustrate that in the domain of 
negative k‘^ (between two curves of Fig.vortices (better 
to say regions with closed current lines) occupy larger 
areas at the same field than in the case k^ > 0 of Fig.[^ 
Hence, vortices appear at the bottom of the suppressed 
order parameter valley. They have normal cores in a 
sense that A = 0 at the center of each vortex and the 
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phase changes by 27r if one circles the center. Still, they 
differ from their Abrikosov “brethren”. The order pa¬ 
rameter changes differently with the distance from the 
center along x oi y directions. Unlike Abrikosov’s case, 
one cannot define the core size as the distance from the 
center to a place with depairing current density. Besides, 
the self-energy of these vortices should be quite small 
because they appear in the region where the order pa¬ 
rameter is strongly suppressed even in zero field. The 
magnetic field is practically constant in the N layer (it is 
exactly constant within our model). In particular, this 
implies that methods of observing vortices based on de¬ 
tecting the vortex field (decoration or scanning SQIUD 
microscopy) will probably not work. On the other hand, 
STM that probes the order parameter value should dis¬ 
cern zeros of A. In fact, the recent STM data show vor¬ 
tices between superconducting Pb islands separated by 
the normal wetting layer [ 1 ]. 

There are many questions remain on properties of vor¬ 
tices within domains of proximity induced superconduc¬ 
tivity. Currents through the SNS sandwich in magnetic 
field should cause vortex motion. Is this motion over¬ 
damped and if it is, what is the drag coefficient? The 
Bardeen-Stephen formula is unlikely to work since it is 
not even clear what plays the role of the vortex core size 
in the normal metal. 

An interesting question concerns superconducting 
fluctuations in the N phase. According to pioneering 
results of Schmid [9] and Prange uni based on linearized 
GL equation, the diamagnetic susceptibility Xd Ihe 
normal phase is proportional to In particular, in zero 
field, Xd diverges as T approaches Tc from above. Here, 
a method to evaluate is offered for any place 

in the {H, T) phase diagram. It would be of interest to 
look at possible differences in diamagnetic susceptibility 
within the region where = — 1 /C^ is negative (between 
the curves of Fig.Q and out of it where k‘^ is positive. 

The author is grateful to V. Dobrovitski, L. Bulaevskii, 
S. Bud’ko, R. Prozorov, P. Canfield, D. Finnemore, 
M. Hupalo for helpful discussions. The Ames Lab¬ 
oratory is supported by the Department of Energy, 
Office of Basic Energy Sciences, Division of Materials 
Sciences and Engineering under Contract No. DE-AC 02 - 
07CH11358. 


Appendix A. Eor the numerical work the integral 
is rewritten to account for the branch point at 
s = l/y/u: 


S = 



(1+^y 

(1 - 


erfc^= — cos(7rcr) erfc 
^/u 


1 

r]\/u_ 

(15) 


Eor the calculation of Hc 2 ,n{T), it is convenient to mea¬ 


sure length in units of hv/27TTcN- Then, we have: 

q 2_ H 

A-f t(2n-h 1) 27rT^^0o Hc2,n{^) 

where A = h/27rTcN^ is the non-magnetic scattering 
parameter, i 7 c 2 ,Ar( 0 ) is the zero-T clean limit upper 
critical field, and t = T/Tcn- 

Appendix B. The assumption of a finite Tcn in the 
main text is in fact not necessary. However, the formal 
treatment of the case Tcn = 0 should take into account 
that A = 0 when the effective coupling is zero. Neverthe¬ 
less, proximity with S results in non-zero Green’s func¬ 
tions F{uj) in the normal metal. This leads to f3 — S = 0 
and to different exponential decay lengths of F for differ¬ 
ent uj = 7 rT( 2 n + 1)/^- The longest length corresponds 
to n = 0, i.e. to u; = ttT, so that calculating the depth 
of pairs penetration one can disregard all n 7 ^ 0 lain]- 



FIG. 7. (Color online) The curve k — 0 for Tcn — 0. Note 
that this curve holds for any mean-free path; the actual tem¬ 
perature and field are T = /ito/27rT and H = c^oho/27 tF. 
Hence at (H^T) plane, on approaching the clean limit, this 
curve shrinks to the origin so that k > 0 everywhere. On 
the other hand, the domain of /c < 0 expands with increasing 
scattering. 


Since in this situation there is no standard energy scale 
related to A or Tc (and no length scale hv/Tc)^ one can 
use the following reduced temperature and field: 


27rr 


T, 


ho = 


(po 


(17) 


In these variables, (3 = 1 + to and u = ho/(l + to)^. 
To find /c(to, ho) one has to solve 1 + to = S{u,a) with 
S taken at n = 0 . Consider, as an example, the curve 
h = 0 at which a = —1/2. Using the form ( p!^ , we have 
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This integral is expressed in terms of generalized hyperge¬ 
ometric functions, which are easily treated with the help 
of Mathematica. Solving numerically 1 + to = So{u,a) 
one obtains the curve of Fig.[^ 
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